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Abstract 

We study the quantum dynamics of diatomic molecule driven by a circularly polarized resonant electric 
field. We look for a quantum effect due to classical chaos appearing due to the overlapping of nonlinear 
resonances associated to the vibrational and rotational motion. We solve the Schrodinger equation associ- 
ated with the wave function expanded in term of proper stationary states, \n)d§\lm) (vibrational®angular 
momentum states). Looking for quantum-classic correspondence, we consider the Liouville dynamics in 
the two dimensional phase space defined by the coherent -like state of vibrational states, and it is found 
some similarities when the quantum dynamics is pictured in terms of number and phase operators. 

PACS: 33.40.+f 

1 Introduction 

The study of quantum dynamics in the interval of parameters where classical chaotic behavior occurs is 
what we call "Quantum Chaos, Chaology, or Quantum Manifestation of Chaos" [li which deals with some 
type of quantum manifestation of the classical chaos, mainly associated with the statistical properties of 
near neighbor levels of energy of the system pQ. In contrast, for quantum systems associated to non chaotic 
classical ones, it is mostly believed that classical dynamical behavior must occur for large quantum numbers 
or high value of the action variable [2] (Rydberg states). In particular, studies of dynamical chaos in atomic 
and molecular systems has been of great theoretical and practical interest [3]- [12] since not enough integrals 
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of motion are found either in the classical or in its quantum system. Different approaches and studies have 
been used for the classical [T3J-Q1)] and quantum (quasi-classical region) [I|)]-[TB] cases, and most of them are 
based on the Morse potential as the inter-atomic interaction [TJj]. On the other hand, the classical study of 
the dynamics of atomic and molecular systems has shown that, under certain conditions, these systems are 
capable of exhibiting a chaotic behavior, even in the case the system has few degrees of freedom. In what 
follows of this introduction and to have a better perspective of the problem, we summarize what Berman et 
al [20] did for the classical part of the problem. 

It is known that the dynamics of a diatomic molecule in a resonant circularly polarized electric field can 
be modeled by the Hamiltonian 

H(I, #;p e , 6; Pv ,<p) = H^ lb) (1, 0) + H<f ot) fa,0;p v , tp) + H mt (I, 0, d, t) (1.1) 

where Hq^ describes the vibration of the molecule along its axis in terms of the action variable / and 
its conjugated angle variable 0, Hq 1 " ^ describes the rotation of the molecule around its transversal axis of 
symmetry in terms of spherical coordinates (r,0,(j>), and Hi nt = —d ■ E(t) describes the interaction of the 
electric dipole moment (d) of the molecule with and external electric field (E(t)). These term are given 
explicitly by 

H? ib \l,0) = ujJ-hx^J 2 , 4 rat) (peAp v ,<p) = A (po + > (1-2) 

z/j,r z V sin I 

and 

H int = -E d - E ° eeff , / _ sin 6 [cos($ - ip + tot) + cos(t9 + ip - uit)] . (1.3) 
2 y /j,uj e 

For the derivation of these expressions, the motion of the molecule has been done with respect the center of 
mass (R = (mxfl + m^f^) j '(mi + m 2 )) and relative ( f = — f\) coordinates associated to the diatomic 
molecule, where \i is its reduced mass. The parameter x e is defined as x e = (h/2u e )(d 2 H ( mb) /dI 2 ). The 
electric field has been chosen as E(t) — Eo(cosu>t, sinwt, 0), and the magnitude of the electric dipole moment 
has been given by d — do + e e f t y/2I / [iio e cosi9, where do = e e //T"o, being e e tf the effective charge of the 
molecule and r represents the point of the minimum on the Morse potential [19] which simulates the atom- 
atom interaction in the diatomic molecule and has been taken up to second order. The average small vibration 
oscillation around the equilibrium point r a is just J Q 2w (Ar) 2 d$ = 2/ cos 2 $//iuj e , with oj e representing the 
angular frequency of the oscillation of the molecule at first order. The dynamical system described by this 
Hamiltonian close to resonance (uj ~ cj e ), and under the condition uj 3> \4>\, \6\ has the following constant of 
motion 

p v — I = constant = Hk, (1-4) 
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The Hamiltonian (1.1) can be written in a more suitable form through a change of variable defined by the 
generatrix function i^^j f, $j n > k,p; r ) = ^( n + 1/2)0$ + <p — wt/U) + /ifct/p 4- which are given by 

ij) = d + Lp - ut, 6 = 6, tp= ip (1.5) 

I = h(n+ 1/2), p = hp, p v = hk + I, T = Clt, (1.6) 



and the Hamiltonian in this coordinates is written as 



1\ . ( l x2 



H(n, ip; p, 6; k, cp) = h(cj e - w) ( n + - ) - hx e cj e I n + - ) + /3 



1 / 



p 2 + -^5- ( k + n + 



sm 2 ' 



1¥ / 1 

~ TV n+ 2 sin6 * cos ^' ( 1 - 7 ) 

where the variable "n" is assumed to have continuous values, and the following definitions have been made: 
P = h 2 /2[irQ, and W — Eoe e ffW2H/fiu) e . This Hamiltonian depends only in the conjugate variables (n,ip) 
and (p, 9) since tp is an ignorable variable, and therefore, fc is a constant of motion. In this way, the Hamilton 
equations with this variables define the four dimensional classical dynamical system 



dn W 1 . „ . . 

dtp . / 1\ 2(3 ( 1\ W . n 

— = Muip — uj) — 2nx„uj f , [ n H — M o — [ k + n H — , sm 9 cos w, 

dr e \ 2 J sm 2 9\ 2 J + 1/2 

2/3 ^, 1\ 2 / 1 n 



rfr sin 3 6» V'" ' 2 J """" ' 2 V ' 2 
f = 20p. (1.8) 

This system has its critical points at p = 0, 9 = it /2, ip — mix with ra£Z, and n = rij with i = 1, 2, 3 the 
roots of a third order polynomial. In the example given by the reference |20j . the parameters associated to 
the diatomic molecule GeO [IT are used, 

huj e = 985.8 cm -1 , H((v e — w) = 15 cm -1 , Hx e oj e = 2.2 cm -1 , /3 = 0.48 cm -1 , 

d = 3.28D, r = 1.621, fj, = 13.1 amu, (1.9) 

where units have been selected such that Ml = IfcfeT = 1 cm^ 1 and T = 0.956 K, and note that they corre- 
spond to a close resonance. Then, there are one center points at ip = tt and n\ rs 2 which forms one of the 
first resonances of the system (for W — 0.05 cm~ x and k = 0), and other resonance which is located near 
n% ~ 1.5 is due to 9 degree of freedom rather that a critical point of the system. The resonances overlapping 
Chirikov's criterion [ST] for appearing of chaos was verified at W = 0.177 cm^ 1 , and total chaotic behavior 
is observed after W = 1.03 cm -1 , see Figure 1. This result points out that classical chaotic behavior appears 
within the first two exited states of the associated quantum system. However, the correspondence principle 
[5] tells us that the quasi-classical behavior of a quantum system is gotten for very large quantum numbers. 
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Therefore, one would not expect any quasi-classical behavior for ground and first exited states of the quantum 
system. 



(a) (b) 




Figure f Poincare map for 8 = ir/2 with 8 > 0, k = 0, 8 = l,and po — and for: (a) W = 0.048 cm 1 ; 
(b) W = 0.68 cm' 1 ; (c) W = 1.03 cm" 1 . 

In this way, in this paper we study the quantum behavior of this system in the region of parameters where 
this classical chaotic dynamics appears. This behavior could be important in the study of diatomic molecular 
clouds for the star born formation or supernova wind shock studies from dying stars |26j , [27] . We solve the 
associated Schrodinger equation, assuming the wave function is a linear combination of the stationary states 
with time depending coefficients and solving numerically the resulting equations for these coefficients , and 
picture the expectation values of the quantum variables in a phase space-like to look for a similarity with the 
classical behavior. 
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2 Quantum dynamics 



2.1 Quantum Hamiltonian 



Our goal is to solve the Schrodingcr equation, 



where H is the Hermitian operator associated to the Hamiltonian, 

1 ( r 



p v 



2/jr 2 \™ sm9 2 I 
For this propose, the following operators are assigned to the observables, 



Eoe e ff 21 . 

— * / sin cosfw + (f — Lot) 

2 y )ILU e 



(2.1) 



(2.2) 



h ( aa^ + 



sm 



p — ► ^ (2.3) 



where a and = (ay represent the usual ascend and descend operators of the quantum harmonic oscillator. 
If {\n}} and {|^m)} are the basis for the harmonic oscillator and the angular moment operators such that 



a\n)=^/n\n), a'\n) = *Jn\n + 1), <va\n) = fi\n) = n\n), [a, a'] = l 



and 



L 2 \lm) = h 2 l{l + l)\lm), L z \lm) = hm\lm), -I < m < I 
the action of d is defined in terms of the phase operator [2"2"] 



(2.4) 



(2.5) 



as 



1 



= a, e 



« . = ( e -«)t = a t 



1 



sun; := |[e«- e"*" 



These operators have the following properties, 



n, cos 



i sin?9, 



n, sin 



-icosd , e lS \n) = \n+ 1), e"^|n) = |n - 1). 



(2.6) 



(2.7) 



From Eq. (2.2 1 and the definitions (2.3 1 - ( 2.5 1, we get the quantum Hamiltonian of the form H = Hq + W, 



where Hq and W are given by 



H = huj e (^a) a + — hx e 



uj„ a 1 a + 



P L 2 



(2. 



and 



W 



W 



cos $ cos((p — Lot) — sin sin(<£ — cot) 



cos i? cos(<£> — out) — sin sm(op — uit) 



(2.9) 



To construct the Hcrmitian operator W, we have used the fact that for any operator A, the operator (A+A')/2 



is Hermitian. Using the commutation relations of Eq. (2.7), one gets the commutativity of the vibrational 
and rotational operators , [/(a, aJ), g{9, tp)] — 0, and we see that hk = L z —h(h + l/2) is a constant of motion, 
that is 



[H, hk} 



II . L - M n + ~ 







(2.10) 



which is the quantum analogue of Eq. ( |1.4[ ) and has the physical meaning that external electric field's photons 
excite the rotational and vibrational degree of freedom with the same number of quanta. The main reason 



for choosing the phase operator as Eq. (2.6) was to be able to get this quantum constant of motion correctly. 



2.2 Time evolution equations 

Using the number states \n) and the spherical harmonic states \lm), we see that Hq is diagonal in the basis 
{\nlm) = \n) ® \lm)}, and the eigenvalues are given by 
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H \nlm) = E n i\nlm), E nl = huj e ( n + ^ j - Hx e oj e (n + - ) + f3l{l + 1) 



(2.11) 



where there is a (21 + 1) degeneration due to quantum number "m". On the other hand, since the relation 
E n +i,i — E n i > must be satisfied, the quantum vibrational number "n" is bounded [19j in the following way 



< n < 



1 1 



(2.12) 



where [x] means the integer part of the real number "x". Therefore, the spectrum is finite. Let us propose 
the solution of Eq. (2.1) of the form 



-itE„,i/h 



\nlm) 



(2.13) 



Now, substituting this equation in (2.1 ) and using the orthogonality relation (n'l'm'\nlm) — 5 n ' n 8i>i5 m ' rn , we 
get the system of equations for the coefficients as 



ihi) ,,, , - V e it( - E "-'''- E "- l)/h D , W ,„ , , 

li v±s n > [' m > t_ ±j nlm ' y n I m ;ntm 3 

nlm 

where the matrix elements W n >i' m > \ n im — {n f l f m'\W\nlrn) are given by 

^Wn'l'm' .nlm 



(2.14) 



W 

Y 



n'\ \ I h + - cos S\n) + (n'\ cos -dJh + - \n) 



(;/' . h + - sini?|n) + (n'\ sin-^y ft + -\n) 



(I m | sin0cos(<^ — ujt)\lm) 
{I'm! | sin 8sijx{(p — cot) | Im) > 



(2.15) 
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From the expressions Al to A7 of the appendix, this matrix elements can be written as 



W n , Vm ,, nlm = - ~ (Vn + 1/2 + y/n' + 1/2) - 5\7 l ) x 



x ( e -^ts2+ 1 s rn+1 — — 



n m 



c lm (2l + 1) 



n' m' 



c Vml {2l' + l) 



(2.16) 



Substituting this expression in Eq. (2.14) and using the same dimensionless variables defined in the intro- 



duction, we get the time evolution equation of the coefficients as 

W 



,-n' _ _ 



(Vn + 1/2 + ^/n - 1/2) — "z!r~7r fl(i-l,l-l,m- J 



c /m e lTf W-),(+) 



\ ' Cl+l ,m-l 



(21 + 3) 



-On-lJ+l,m-l 



+ 1/2 + V" + 3/2 /0 , , n Ai+1 ,m+l 



c /m (2/ + l) 

«"fi n! , (+),(-) 



+!,;+!, m+l 



(2.17) 



where we have made the definitions D 1 = dD/dr, fi n i,(±),(-) = E n ±\^\ — E n i =F hu, and Q n i,(±).(+) 



E n ±ij+i — E n i =F Huj. The selection rules deduced from (2.17) are 



Al = ±1, An = ±1 and Am = ±1 



(2.18) 



Note that the last two terms of these expressions are a consequence of the constant of motion (2.10). The 



time evolution of the coefficients in Eq. (2.17) and the selection rules in Eq. (2.18) are similar to the electric 



dipolc transitions in an atom, except with the extra selection rule of n. Furthermore, suppose we are initially 
in a given state 1^(0)) = |no^o m o) & n d we set the frequency u to be such that it is almost in resonance 
with the frequency of an allowed transition, say \nflfmf) (that is Hu) ~ E nf j f — E no j ). For this case and 



neglecting the non- resonant transitions, the equations of motion (2.17) becomes 



iD' = ae irUr D f and iD' f 



ae -Do, 



(2.19) 



where a and Q r are defined as a = — (^y/nQ + 1/2 + y/nj + 1 J2j Wci 0imo /Aci ftmf (2lf + 1) and Q r = E nft i f 
£ n 1 — huj. In matrix notation, Eq. (2.19) is written as 



dr 



D, 
D 



a 







if) 




which in terms of the Pauli operators becomes 



i — = a[cosf2 r T a x +sinf2 r T a y ] with \tp) = 

dr 



D f 
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and this one is of the form 
. d 



= H at\4>), with H at = a 



cos O, 



r er T + sin Cl r T a 



(2.20) 



which is the Schrodingcr equation for a two level atom introduced in a circularly polarized electromagnetic 
field 0. 



2.3 Numerical results 



We solve numerically the Eq. (2.17), considering only the coefficients D n i m for n, I < 3. We use the same 



parameters used in the classical numerical calculations, Eq. (1.9), which implies to have a close resonant 
transition between the states 1 100) and |211) with |Oxo,(+),(+) | = 0.82. Higher order of excitation are not 
considered since we want to see what happen the states closely related with the classical ones, where classical 
chaos appears. The results of the numerical simulations are shown in Fig. [2j» and Fig. |}» . The Fig. 
[2^ a) shows the transition probabilities for small values of W before there is a considerable mixing of states 
(observe that |c2ii| 2 < 0.5). Fig[2^b) shows the same probabilities but for the value W — 1.03 cmT 1 , and 
here one can observe a considerable overlapping of the main transition coefficients involved in the dynamics 
(this case corresponds to have classical chaotic behavior). 



f \ W = 0,048 

cL I --- W = 0.19c 
K ' W - 0.68 c 



(b) — W 





v< W W \.' i W i \, \/ V/ 



15 20 25 30 35 40 45 50 



5 10 15 20 25 30 35 40 45 50 




1.03 
1.02 




5 10 15 20 25 30 35 40 45 50 



10 15 20 25 30 35 40 45 50 



Figure 2 Time evolution of the total probability, the probability amplitude of the state 1 100) (upper 
line) and the state |211) (lower line) for different values of the perturbation: (a) W — 0.048 cm -1 , solid; 
W = 0.19 cm- 1 , dashed; W = 0.68 cm" 1 , bold; (b) W = 1.03 cm -1 . 



Note that the value of W for classical chaos to appear is W c h = 0.177 cm -1 , and this one does not coincide 
with the value for mixing of states (W « 0.68 cm -1 ), although they are of the same order of magnitude. 
We see also that the classical value of the closed classical resonance suggests overlapping between quantum 
states in n = 1 and n = 2, as we precisely observe in our simulations, which is consequence of the resonant 
transition frequency between the states 1 100) and 1 2 1 1) . 
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2.4 Quantum phase space pictures 



In this section we try two different approaches to get a better relation between the quantum and classical 
dynamics. The first and most used approach [22j- L 24J, is to use the phase space representation in terms of 
the expectation value of the dimensionless canonical variables X and P 



X = 



a + a) 



P = 



(2.21) 



2 ' 2i 

The results of the numerical simulations of this approaches is presented in Fig. [3] , where we used the same 
parameters of expression (1.15) and W = 1.03 cm -1 . For the initial state wave function of the system we 
chose a poisson-like distribution in the coefficients D n i with the maximum in value in the state -Dioo- This 
initial state is determined by the following coefficients 



^ooo = 
D ioo = 
-D200 = 



'1.5 
12' 

IT 
12' 
T5 
T2 



I>110 = 

-D210 = 



/0.2 
12' 

fOA 
12' 

/02 
12' 



-D020 = 
-D120 = 
-D220 = 



0.05 

IT 

oT 

12 

O05 

~Y1' 



(2.22) 



Based in the properties of the coherent states of the harmonic oscillator, this selection not only gives a well 
defined initial value of the expectation values, but also permits a further study in terms of the Liouville 
dynamics for both the classical and quantum case |25j . The initially values with quantum number m = is 
allowed by the slow oscillations of 9 and ip assumed in the model. Although the dynamics of each variable 
seem to be stable and similar to each other, the phase space representation does not seem to give any pic- 
ture alike to the classical dynamics of the system (see reference [20]), i.e., the phase space in terms of the 
canonical variables (X) and (P) shows no resemblance to the classical dynamics, neither it seems to suggest 
any transition to quasi-classical chaos. 
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In the second approach, we will use the phase space representation in terms of the expectation value of the 
phase $ and the number operator n. In polar coordinates, the (h) will correspond to the radius and arg ( e% j 
to the angle of the phase space representation of this set of variables, h — a 1 a and e vd = (1/vaaJ). The 
expectation values of these variables represents the classical analogous of the variables on Figure 1 above. 
In the upper left plot of Fig. [4] it is shown the time evolution of (n) which resemblances to the classical 
case in terms of the main and different frequencies with which it oscillates. For the numerical simulations 
results presented in these figures we used the same parameters as in the section [2^31 and the same initial state 



(2.22). The phase space picture in terms of the operators n and e seems to have a little bit resemblance 



with the classical results, perhaps because the dynamics of n resembles the classical part. Also, the sudden 
slow changes of arg i^ 1 ^ seem to suggest some kind of relation with the resonances of the classical case. 




Figure 4 Phase space like picture with the same parameters as in the classical case in terms of the 
expectation values (n) and (arg(e 1 ' 6 )). 
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3 Conclusions and comments 



We have studied the quantization of a diatomic molecule by solving the Schrodinger equation with the 
known Hamiltonian of the diatomic molecule with a circularly polarized resonant rf-field, written in spherical 
coordinates (rotations) and angle-action variables (vibrations). The wave function was expanded in a finite 
combination of a proper stationary basis with time dependent coefficients, and the system of equations for 
these coefficients was obtained. Using the same parameters as in the classical case, a near resonant transition 
between the states 1 100) and 1 21 1) is gotten, which correspond to the closer integer numbers for n where the 
classical non linear resonances appeared, n\ ~ 2 and n 2 ~ 1.5. The critical value for the mixing between this 
two states to occur is within the same order of magnitude as the classical value for resonances to overlap. Using 
as an initial state a poisson-like distributed wave function in the quantum numbers n and 1 with maximum in 
the resonant state 1 100) , we try two different approaches to see the quantum phase space expectation value 
dynamics and compare it with the classical case. The usual approach, using the canonical variables X and 
P, fails to provide any intuitive picture of the classical case. On the other hand, the approach using the 
expectation value of e 1 ^ and n suggests some resemblance and relationship with the classical case. 

Therefore, we have here the following situation, on one hand, the correspondence principle tells us that 
we must have the quasi-classical behavior for this quantum system at very large quantum number. However, 
due to the quantum dynamics involves its first few states, quasi-classical chaotic behavior can not be obtained 
here, but classical chaotic behavior is obtained just for these first (continues) states. We do not see how this 
quasi-classical limit could happen for this quantum system with double nonlinear resonance. So, as one could 
expect for this case, quantum dynamics does not follow the classical one. 
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5 Appendix 



Matrix elements of some operators 



(n'\\l cos ^l n ) 



n ' + 2 ( S n',n-1 + S n ',n+l) , 



1, , 1 



(n'\cos$\ h+ -\n) = -Jn+- (5 n /, n _i + S n > >n+1 ) , 



2t 



[u'\ sin ^\/n + 



(n'U/n+ -sini?|n) = -i/n' + - (£„',„+i - 5„', n -i) 



2/ 



n + 2 (^"',"+1 — <^n',n-l) : 



(Z'm'| sin#cos<£|Zm) = 

(Z'm'l sin^siin£|Zm) = 
The cocfFicients a m are given by 



cm+l cm— 1 



rm+1 

21 + 1 C; 



+ 1 C V . 
rm-l 

, m [dl ' > 21' + 1 c Vm , V 1 



Clr 



2(1 



(2Z + l)(Z-ra)! 
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